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Abstract. Using a Tsallis nonextensive approach, we analyse distributions of transverse spectra of jets. 
We discuss the possible influence of energy conservation laws on these distributions. Transverse spectra of 
jets exhibit a power-law behavior of l/Pr with the power indices n similar to those for transverse spectra 
of hadrons. 
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1 Introduction 

For some time now it is known that transverse momentum 
spectra of different kinds measured in multiparticle pro¬ 
duction processes can be described by a quasi power-law 
formula Q 

hiPT)=c{i+^y\ ( 1 ) 

which for large values of transverse momenta, px ^ nT^ 
becomes scale free (independent of T) power distribution, 
1/py. This was first proposed in [TJ[5] as the simplest for¬ 
mula extrapolating exponential behavior observed for low 
Pt to power behavior at large px- At present it is known 
as the QCD-inspired Hagedorn formula [51H] . Distribu¬ 
tion o, nowadays, has been successfully used for a de¬ 
scription of multiparticle production processes in a wide 
range of incident energy (from a few GeV up to a few of 
TeV) and in broad range of transverse momenta (see, for 
example, reviews mm)- In particular, it turned out that 
it successfully describes a transverse momenta of charged 
particles measured by the LHC experiments, the flux of 
which changes by over 14 orders of magnitude [iniiiniin]. 

In many branches of physics Eq. o , with n replaced by 
n = 1/ (g— 1), it is widely known as the Tsallis formula [121 
nm]. In this form, Eq. m is usually supposed to rep¬ 
resent a nonextensive generalization of the Boltzmann- 
Gibbs exponential distribution, exp{—px/T), used in a 
statistical description of multiparticle production processes 
Such an approach is known as nonextensive statistics m 
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^ Throughout the whole article a natural units (Boltzmann 
constant and speed of light) are used by setting its values to 
unity. 


mm in which the nonextensivity parameter q summarily 
describes all features causing a departure from the usual 
Boltzmann-Gibbs statistics. In particular it can be shown 
that the nonextensivity parameter q is directly related to 
the fluctuations of the parameter T identified with the 
’’temperature” of the hadronizing fireball [TCIITB] . 

It was shown there that such a situation can occur 
when the heat bath is not homogeneous and must be 
described by a local temperature, T', fluctuating from 
point to point around some equilibrium value, T. Assum¬ 
ing some simple diffusion picture as being responsible for 
equalization of this temperature [5l ll5l[T6] one obtains the 
evolution of T' in the form of a Langevin stochastic equa¬ 
tion with the distribution of I/T', g(l/T'), emerging as 
a solution of the corresponding Fokker-Planck equation. 
It turns out that in this case g(\/T') takes the form of a 
gamma distribution. 


ff(i/r') 


nT 

T^) 



n—1 

exp 



( 2 ) 


After convolution of the usual Boltzmann-Gibbs exponen¬ 
tial factor ex-p{—px/T') with Eq. ([2]) one immediately ob¬ 
tains a Tsallis distribution, Eq. o, with 


< r >2 

Var{T') ’ 


(3) 


directly connected to the variance of T'. This idea was fur¬ 
ther developed in [T^f^Oll^ (where problems connected 
with the notion of temperature in such cases were ad¬ 
dressed) . The above makes a basis for the so-called super¬ 
statistics mm- 

Applications of Tsallis distributions to multiparticle 
production processes are now numerous. To those quoted 
previously in [5] one should add some new results from El 
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HHH and also presented in m- The most recent applica¬ 
tions of this approach come from the STAR and PHENIX 
Collaborations at RHIC and from CMS P51[?7] . 

ALICE |281[21I] and ATLAS [20] Collaborations at LHC 
(see also a recent compilations [21112^ 1. 

It is now empirically well-documented that transverse 
spectra of both hadrons and jets exhibit a power-law be¬ 
havior of 1 /p^ at hight pt- This observation (usually inter¬ 
preted in terms of non-extensive Tsallis statistics) meets 
difficulty. Whereas in it has been advocated that 

the power indices n for hadrons are systematically greater 
than those for jets, in [22] the values of corresponding 
power indices were found to be similar, strongly indicat¬ 
ing the existence of a common mechanism behind all these 
processes. 

In Section [2] the influence of constrains, forcing the use 
of conditional probabilities for Tsallis distribution (and 
the case of Boltzmann statistics in Appendix) is discussed. 
In Section [2] we report the results concerning the conse¬ 
quences of the use of conditional Tsallis distribution for 
description of transverse spectra of jets. Section U] is our 
summary and discussion of results. 


2 Conditional probability - an influence of 
conservation laws 



Fig. 1. (Color online) Conditional probability distribution, 
h{Ei\U), as a function of of Ei/U for Tsallis statistics (n = 7 
and T = 1 GeV). 


Let be a set of v independent identically dis¬ 

tributed random variables described by some parameter 
T. For independent energies, each distributed 

according to the simple Boltzmann distribution: 


gi {Ei) = - exp 




Fig. 2. (Color online) Ratio of the conditional distribution 
function h[Ei\U) and a single particle distribution hi{Ei) as a 
function of Ei/U, for Tsallis statistics (n = 7 and T = 1 GeV). 


the sum U = Ei is then distributed according to the 

distribution, gv{U) = [rr{v)]~^{U/TY~^ eyip>{—U/T). 

If the available energy U is limited, for example if [/ = 
Ei = const, then we have the following conditional 
probability for the single particle distribution, g (Ei): 


g{Ei\u) 


gi iEi)gu-i {U - Ej) 
g.{u) 


(5) 


In the above consideration the value of the T parameter 
does not fluctuate. More details is given in Appendix and 
review m- 

Now let us consider a situation in which the parameter 
T in the joint probability distribution 

V 

g{{Ei,...,.}) = X{gi{Efl ( 6 ) 

i=l 

fluctuates according to a gamma distribution, Eq. ([2|). In 
this case we have a single particle Tsallis distribution 


hi {Ei) 


n — 1 
nT 



(7) 


and a distribution oi U = Ei is given by (cf. [22]1: 


K{U) = 


r{v + n-l) fu 
ur{v)r{n-i) \ t 


1-b 




1 — v—n 


.( 8 ) 


If the energy is limited, i.e., if C/ = Ei = const, we 

have the following conditional probability: 


hiEi\U) = 


hi (Ej) {U - Ej) 

hAU) 


(4) 
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X 


(zz- l)(n- 1) {nT + U) ^ 

(n - 2 + 1 /) nTU \ U ) ^ 




E, 


nT + U 


2 — i'—n 


(9) 


For n —)> c» Eq. reduces to Eq. ©• On the other 
hand, for large energy iJJ —>■ oo) and large number of 
degrees of freedom (y —^ oo), the conditional probability 
distribution reduces to the single particle distribution 
given by Eq. ©• 

For Ei U the conditional probability ® can be 
rewritten as 


hmu) 


(zz + 2)(u-l) / 
nT{n — 2 + iy) \ nT J 


( 10 ) 



which, when additionally zz ^ 1, reduces to Eq. ®. 

The results presented here are summarized in Figs. [T] 
and [H which show how large are differences between the 
conditionalTsallis distribution ^(^117) and the usualhi{Ei) 


Fig. 4. (Color online) Jet differential yields as a function of 
pr/y/s or Et/\/s {pr and Et distributions are shown by full 
and open symbols, respectively) fitted by the conditional Tsal- 
'lis distribution given by Eq. 


3 Imprints of constrains in transverse spectra 



Fig. 3. (Color online) Transverse spectra for jets {pr and Et 
distribution are shown by full and open symbols, respectively) 
fitted by the conditional Tsallis distribution given by Eq. ®. 



Fig. 5. (Color online) Jet differential yields as a function of 
Pt/U or Et/U {pt and Et distribution are shown by full and 
open symbols, respectively) fitted by the conditional Tsallis 
distribution given by Eq. 


Now we concentrate on the experimental distributions 
of transverse momenta 


/(tt) 


(Pa 

PTdpTdri 


( 11 ) 


of jets observed at midrapidity and in a narrow jet cone 
defined by i? = y /Ap + A(f>^ (where Acj) and Ap are, 
respectively, the azimuthal angle and the pseudorapidity 
of hadrons relative to that of the jet). Data on p+p and p+ 
p interactions, covering a wide energy range from 0.54 TeV 


up to 7 TeV, obtained by CMS [50], CDF II [57], DO [55] 
and UA2 [55] experiments was used for the analysis (cf. 
Table[T]for more details). Because data, which we analyze, 
are presented at midrapidity (i.e., for ry ^ 0) and for large 
transverse momenta, pp p, the energy E of a produced 
jets is roughly equal to its transverse momentum pp (and 
transverse energy Ep), which we shall use in what follows. 

Transverse spectra for jets f{pp), fitted by the condi¬ 
tional Tsallis distribution h{pp\U) given by Eq. ([9]), are 
shown in Fig. (5] Parameters of the Tsallis distribution 
used to fit the data are shown in Table [5] An available 
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Table 1. Jet production in p + p and p + p interactions. 


Exp. 

Reaction 

Variable 

^ [TeV] 

R 

\v\ 

CMS 

P + P 

PT 

7.0 

0.5 

< 0.5 

CDF II 

P + P 

Pt 

1.96 

0.7 

< 0.1 

DO 

P + P 

Et 

1.8 

0.7 

< 0.5 

DO 

P + P 

Et 

0.63 

0.7 

< 0.5 

UA2 

P + P 

PT 

0.54 

- 

< 0.85 


Table 2. Parameters of the Tsallis distribution used to fit the 
data. 


Exp. 

[TeV] 

n 

m 

T [GeV] 

U [GeV] 

CMS 

7.0 

6.95 

1900 

1.08 

2200 

CDF II 

1.96 

7.55 

1340 

1.08 

1190 

DO 

1.8 

7.7 

1050 

1.08 

960 

DO 

0.63 

8.0 

800 

1.08 

500 

UA2 

0.54 

8.3 

800 

1.08 

490 


energy per degree of freedom 

ks0.22.1„(|5) (12) 

increase logarithmically with energy y/s, while vj = 

12.7 - 2.73 In^ + 0.15 {In^f and U/^ = 2.35 - 
0.23 Iriy/s. The slope parameters decrease with energy as 

n= 12.25 (13) 

and numerically are comparable with n{s) = l/{q{s) — 1) 
for q{s) evaluated for transverse spectra of charged parti¬ 
cles (for charge particles n = 14.78(-ys)“° °®) [30]. 

Of course many other parametrisations, than this given 
by Eq. m, is possible in our limited interval of energy 
(0.5 TeV - 7 TeV). For example n = 4 + 10(-ys)“°'^^^ for 
jets and n = 4-|- 15.6(-yi)“°-^ for charged particles. How¬ 
ever all parametrisations lead to conclusions that rijet = 
nchar at TeV energies (1.2 TeV - 1.4 TeV). 

The flux at given pp (or Et) increase with interaction 
energy ^/s as 

f{pT = 100 GeV) = 5.5 • 10-®(yi)i ®® - 0.6. (14) 

However a comparison of fluxes at transverse momenta 
close to the interaction energy (cf. Fig. |4| show decrease 
with energy 

/(pt/\/s = 0.1) = 10^‘‘(^)“^-'^. (15) 

With increasing interaction energy, jets are produced with 
lower relative transverse momenta and paradoxically for 
lower -y/s we are closer to the kinematical limits. For trans¬ 
verse momenta close to available energy U, presented in 
Fig. [S] fluxes of jets show rapid decrease with pr/U —>■ 1. 

4 Summary 

Transverse momentum spectra of jets, spanning over 12 
orders of magnitude in the observed cross sections, can be 


successfully described by the conditional Tsallis distribu¬ 
tion. Those spectra exhibit a power-law bahavior of 1 /pp 
for which n~7 — 8(nRi8at 800 GeV and decrease with 
energy, n = 8(-v/s/800)“° °®^. Observed deviation from the 
power-law at tail of distribution is caused by kinematical 
constrains. Fluxes of jets for pp/i/s = const are higher 
for lower incident energies, and are comparable for high 
Pt/U where fluxes have been lowering seemingly. 

Comparison of power indices n for jets with those for 
charged particles produced in minimum bias collisions is 
shown in Fig. |B| The values of the corresponding power 
indices are similar, indicating the existence of a common 
mechanism behind all these processes [33] . The suggestion 
that power indices for jets are almost 2 times smaller than 
those for hadrons piTT)] . comes mainly from a different 
parametrisation of hadrons and jets spectra (not a simple 
Tsallis distribution) and can not be accepted by looking 
directly at the corresponding distributions. 


^ i 




d 

o 


- 

[ 1 

jets: 

■ Pt 
• Et 

charged particles: 

' □ CMS 
o WWCT 

0 

s° 

9 

0^ 

10® 


10"^ 


^/7 [GeV] 

Fig. 6. (Color online) Slope parameters n extracted from the 
fits to the transverse spectra of jets (full symbols) in compari¬ 
son with slope parameters for distributions of charged particles 
(open symbols: CMS PCT?7| and WWCT [TT]'l. 


Appendix A: Conditional Boltzmann distribu¬ 
tion 

For independent energies, each distributed ac¬ 

cording to the simple Boltzmann distribution, Eq. (jT]) the 
sum 

U = Y,E, (A-1) 

i=l 

is distributed accordingly to the following gamma distri¬ 
bution, 

(?) (-?) ■ 

If the available energy U is limited, for example if 17 = 
const, then we have the following conditional probability 


































M. Rybczyriski and Z. Wlodarczyk: Transverse momentum of jets ... 


5 


for the single particle distribution (cf. Eq. (O): 


9{E,\U) 


U \ U) 


(A-3) 


This is nothing else than the well-known Tsallis distribu¬ 
tion 


9{E^\U) 



l-(l-g) 


(A-4) 


with nonextensivity parameter 


q = 


V — S 
v-2 


(A-5) 


which is always less than unity and T = U/{v — 1). 

Actually, such distributions emerge directly from the 
calculus of probability for situation known as induced par¬ 
tition [41j . In short: jz — 1 randomly chosen independent 
points {C/i,..., t/jy-i} breaks segment (0, U) into iz parts, 
length of which is distributed according to Eq. (IA-31) . The 
length of the such part corresponds to the value of 
energy Ek = Uk+i — Uk (for ordered Uk)- One could think 
of some analogy in physics to the case of random breaks 
of string in iz — 1 points in the energy space. 

We end this part by a reminder of how Tsallis distri¬ 
bution with q < 1 arises from statistical physics considera¬ 
tions. Consider an isolated system with energy U = const 
and with n degrees of freedom. Choose a single degree of 
freedom with energy E (i.e. the remaining or reservoir en¬ 
ergy is Er = U — E). If this degree of freedom is in a 
single, well defined state, then the number of states in the 
whole system equals 17(1/ — E). The probability that the 
energy of the chosen degree of freedom equals to E is then 
P{E) oc 17(1/ — E). Expanding (slowly varying) In 17(A) 
around U, 


lnl7(C/-A) 


CX3 


E 


1 In 17 

E. dE^'"'’ 


(A-6) 


and (because E <^U) keeping only the two first terms one 
gets InA(A) oc In 17(A) oc —/3A, i.e., P{E) is a Boltzmann 
distribution given by Eq. o with 

_ de/ gin 17 (A^) 

T ^ dEr ' ^ ~ ^ 

On the other hand, because one usually expects that 
{Er) oc (Ar/^)“^‘^~“^ (where 01,2 are of the order of 
unity and we put oi = 1 and, to account for diminishing 
the number of states in the reservoir by one, 02 = 2) [42l 
03] , one can write 


gAfc 


oc (-1)'=A:! 


iy-2 


{-l)^k\ 


(zz — 2)^ 


(A-8) 


and write the full series for probability of choosing energy 
A: 


, , n{u-E) 

m) 


(-1)'= (-/3E)^+^ 

k -I-1 (zz — 2)^ 

.fe=o ^ ' 


= C 




(v-2) 


= /3(2-g)[l-(l-9)/3A]T^; 


(A-9) 


(where we have used the equality ln(l-|-a:) = 

1)]). This result, with g = (zz — 3)/(zz — 2) < 1, coincides 
with the results from conditional probability and the in¬ 
duced partition. 
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